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Abstract 

^ In this paper we consider the problem of determining the boundary perturbations of 

■^Uj I an object from far-field electric or acoustic measurements. Assuming that the unknown 

r) . scatterer boundary is a small perturbation of a circle, we develop a linearized relation 

between the far-field data and the shape of the object. This relation is used to find the 

Fourier coefficients of the perturbation of the shape. 
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^ '. 1 Introduction 



The field of inverse shape problems has been an active research area for several decades. 
Several related scalar problems belong to this field: electric and acoustic scattering form 
two large classes. In direct problems one wants to calculate the field outside a given object. 
In two common situations, one knows either the values of the field on the object (the 
f~^ , Dirichlet problem) , or the values of the normal derivative of the field on the boundary (the 

Neumann problem). Inverse shape problems involve reconstructing the object shape from 
measurements of the electric or acoustic field. Differently from Direct problems which are 
usually well posed, inverse problems are ill posed: the solution has an unstable dependence 



H ! on the input data. 



The formulation of the electric scattering problem is based on the quasi-static approx- 
imation and the related Laplace equation for the electric scalar potential. When a perfect 
conductor is exposed to extremely low-frequency electric fields, the problem is equivalent to 
the Dirichlet boundary value problem for the Laplace operator. 

The sound-soft acoustic scattering problem is characterized by the condition that the 
total field vanishes on the boundary of the scatterer. Thus, acoustic scattering is equivalent 
to the Dirichlet boundary value problem for the Helmholtz operator, with the scattered field 
equal to the negative of the known incident field. 

These two problems are frequently solved by methods of potential theory. The single- 
and double-layer potentials relate a charge density on the object boundary to the limiting 
values of the field and its normal derivative. The resulting integral equations are then solved 
in an appropriate function space, a common choice being the Lebesgue space L^ 
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In this paper, assuming that the unknown object boundary is a small perturbation of a 
unit circle, we develop for both electric and acoustic problems a linearized relation between 
the the far-field data and the shape of the scatterer. Under this purpose, we investigate 
the Dirichlet boundary value problem outside the object entering the Dirichlet data as 
parameters and the shape of the object as variables. 

The linearized relation between the far-field data and the object shape is used to find the 
Fourier coefficients of the boundary perturbation of the object. Suppose that the angular 
oscillations in the perturbation are less than 1/n. In order to detect that perturbation, it 
turns out that one needs to use the first n eigenvectors of the Dirichlet-to-Neumann operator 
corresponding to the unperturbed shape as the Dirichlet boundary data. We may think that 
this result is quite general. When the unknown object is a C^-perturbation of a disk, we 
obtain asymptotic formulae for the Dirichlet-to-Neumann operator in terms of the small 
perturbations of the object shape, and it is worth mentioning the expansions of Dirichlet- 
to-Neumann operators for rough non-periodic surfaces f5| HI and for periodic interfaces [S] . 

Our approach relies on asymptotic expansions of the far-field data with respect to the 
perturbations in the boundary, in much the same spirit as the recent work p and the text 
[1]. We consider only the two-dimensional case, the extension to three dimensions being 
obvious. In connection with our work, we should also mention the paper by Kaup and 
Santosa [B] on detecting corrosion from steady-state voltage boundary perturbations and 
the work by Tolmasky and Wiegmann [10' on the reconstruction of small perturbations of 
an interface for the inverse conductivity problem. 

We deal with electric problems in section 2 and 3, and acoustic problems in 4 and 5. 

2 Formulation of the Electric Problem 

We consider the reconstructing problem of the perfect conductor D^ which is the small 
perturbation of the unit disk D described by a Lipschitz function / and a small scale factor 
e, that is 

dD,{=dD + efee):=[{l + ef{B))ee, 0e[O,2^]}, (2.1) 

where ee — (cos 6*, sin 6*). 

2.1 Electric Scattering problem 

If we apply an initial potential v^ to E^ which is homogeneous except for the perfect conduc- 
tor De, then the derived electric potential v is given hy v = v'^ -\- v'^ , where v" is the solution 
to _ 

' Av" = 0, in M^ \ £)^^ 

< w'* = -?;' + C(constant), on dD^, (2.2) 

v'^{x) — > 0, as |a;| -^ ao. 

We denote Vq as the perturbation of electric potential due to the conductor D, i.e., 

Av^ = 0, in R2 \ D, 

-y^ = -w'-f C(constant), on dD, (2.3) 

^0(2;) — > 0, as |a;| -^ cx). 

By obtaining a linearized relation between ef and the e-order term of (u* — Vo)(r, 0) as 
r — > 00, we try to recover D^. 



Au = 0, 


in M^ \ L»e, 


u{l + ef{9),e) = ^{d), 


for 6*6 [0,27r], 


u{x) = constant + 0(1/ a; ), 


as |a;| — > +oo. 



It follows from the Taylor series expansion of v^ near D that 

v\l + ef{e),e) = «^(l, 0) + efdrv\l,e) + 0{e^). (2.4) 

Here we used the polar coordinates x{r, 6) ~ {rcosO, r sin6'). We investigate v'' by consider- 
ing two exterior boundary value problems, one with Dirichlet value (— w'(l, 6)) and the other 
with {~efdrV^{l,9)) on dD^. To do that, we formulate the fixed boundary value problem. 

2.2 Fixed Dirichlet boundary value problem 

When the boundary value is prescribed as the 27r-periodic function 'J, the voltage potential 
outside the conductor D^ is given by the harmonic function u which satisfies the following: 



(2.5) 



We let uo be the voltage potential outside the unit disk D with the fixed Dirichlet data 
^ on the boundary, i.e., 

Auo = 0, in R^ \ D, 

uoil,0) = *(6'), for e e [0, 27r], (2.6) 

uo{x) — constant + 0(l/|a;|), as \x\ —* +oo. 

We obtain the linearized relation between the boundary interface of the conductor D^ 
and {u — uq) at infinity, especially when ^ is given by a C^-function or a Lipschitz function. 

3 Electric far-field formula and Inversion algorithm 

3.1 Linearized relation for the Dirichlet problem 

We start by explaining the main idea to obtain the linearized relation. 

To derive the asymptotic expansion of the solution u to (|2.5p with the given boundary 
data \1/, we apply the field expansion method (F.E) (see ^9j). Firstly, we expand u in powers 
of e, i.e., 

w(r,0) = ^u„(r,0)e". (3.1) 

11=0 

Now expanding in terms of r and evaluating (|3.ip at r = 1 + e/, we obtain that 

u{l + ef{e),e) = uo(l, 0) + e(ui(l, 9) + 9,uo(l, 0)f{e)^ + 0(e2). 

Since u(l + ef{9),6) and ^0(1,6*) have the same value, ui can be considered as the 
decaying harmonic function which satisfies 

u,{l,0) = -Afo{^){0)f{e), (3.2) 

where 

Mo{-i'm^drUoil,9). 



In other words, A/q is the Dirichlet-to-Neumann operator of _D, and it can be expressed as 
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A/'o(*)(0) = - y^ [na„(^) cos7i6> + n6„(^) sinnO 



(3.3) 



where a„('I') and 6„(^) are the fourier coefficients, that is 



1 



2tt 



2Tr 



a„(*) = -/ ■9{9)cos{n9)d9, bni'i/)^- *(e')sin(» 



From (j3.2p and the expansion of the harmonic function outside a disk, we have 



+00 
ui(r,0) = -^— [a„(AAo(*)/)cosn0 + 6„(AAoW/)sinn0 

n=0 



r > 1. 



Therefore, 



(u - uo){r, 6) ai(AAo(*)/) cos0 + Si(AAo(*)/) sin0 



+ C, r>l. 



where C is a constant. More precisely, we have the following theorem and give the proof in 
Subsection [ 



Theorem 3.1 For a 2'K-periodic function ^, we let u and uo he the solution to p.Sp and 
p.6p . respectively. 

1. Let * e C"*([0, 27r]). for r > 1, we /laue 

(M-wo)(^e) = --fai(A/'o(*)f)cos0 + 6i(AAo(*)/)sin0l+C + O(ei/r + e/r2), (3.4) 
r L J 

where C is a constant, and 0(e2 /r + e/r^) depends on the Lipschitz constant of f and 

ll*llc^- 



2. For a Lipschitz function ^ , we have that 

{u - uo){r, 6) = C + O0lr + l/r^), for r > 1, 



(3.5) 



where C is a constant, and 0{e2 /r+ 1/r^) depends on the Lipschitz constant of f and 

Remark 3.2 For the case of C^ -perturbation of the interface, i.e., f G C^([0, 27r]), the 
error term of (j3.4p and (j3.5p can be replaced by 0{e^ /r + e/r^) a^ii 0{e/r + 1/r^). 

In connection with the results for rough non-periodic surfaces |8l |4] and for periodic 
interfaces [5] , we expand the Dirichlet-to-Neumann operator AC/ of D^ which is defined by 



1 unit 
Note that fj, is given by 



du 



(y), y=il + ef{0))e 



where Vy is the outward unit normal vector to Z?c 



Ne 



where 



'~\Ne\' 
Ne = il + efi9))ee - efre, rg = {- sin 9 , cos 9) . 



(3.6) 



Here / is the derivative of / with respect to 6. From the fact that 

1 



\Ne\ 



= l-ef + Oie'), 



(3.7) 



it follows that 

Kfi^m = (1 - ef){S/u,Ne) + 0{e^) 



(l-^/)|(l + ^/)§^ 



ef du 



du 
dr 



xdu 

r=i+efie) do 



r=i+£/ 1 + ef 89 

r=l+e/(e) 



r=l+£/. 



+ 0{e^) 



Applying (|3.ip . we obtain 

Nef(-^m ^ dru^{l, e) + e(drui{l, 9) + dluo{l, 9)f{e) - deuo(l, 9)f{0) 
Defining an operator Pq by 
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Vq{^)(9) := - ^ \{n + l)a„(*) cos 716* + (n + 1)6„(*) sin n6l 



we have 



92^0(1, e)=2?oAAo(*)(0). 



(3.8) 
(3.9) 



Lemma 3.3 For f e C^{[0, 2tt]) and * E C^{[0, 2tt]), 



AC/(*) = AAo(*) + eAA/(*) + 0(ei), 



where 



We give the proof in Subsection [ 



(3.10) 



3.2 Algorithm for the Inverse Shape Problem 

For an entire harmonic function w% we let v'^ and Vq be the solution to (12. 2p and (|2.3p . 
respectively. The Dirichlet values of the solutions are given by 

Vo\dD = —v\l,9) + constant, 
v'loD, = -w'(l, 9) - ecf^wXl, 0)f{O) + constant + O(e^). 



Note that 



drv\l,9)^-Uo{v'\oD)- 



Here we have the minus sign on the right hand side because Nq is the Dirichlet-to-Neumann 
operator for the exterior harmonic functions. Applying (13. 4p and (|3.5p with letting ^ = 
—v^\dD and ^ = A/'o(w*|aD)/, respectively, we obtain for r » 1 that 



iv^-v^„){r,9)^2- 
r 



ai{Afoiv'\aD)f) cos9 + h(Afoiv'\oD)f] sint 

Now define entire harmonic functions w"'* and w"'*, for n G N, by 

i;"'Vr,6l) = --r"sinn6l, w"'Vr,6l) = -r"cosn6l. 
n n 



(3.11) 



Let v^'^ and uj"^^ be the solution to (|2.2p with the initial potential w"'* and w"'% respectively. 

Q and Wq 



In the same way, define Uq'* and w^'^ as the solution to (|2.3p . Let 






e 
From (|3.1ip . it follows that 



(w"'*) - 2ai(f-Afo{--sinne)) ^ 2ai(f ■ smnOY 



Cl 

n 



di{v"'') - 26i(/-A/'o(--sinn6')J = 2&i(/-sinn^ 

By the same way, we obtain 

ci(w"'*)--2ai(/-cosne'), 
di{w"^') ^ -2bJf ■ cosneY 

Thus we obtain that 

1 n^TT 



ci(w"'*)±di(w"^*) == - / 2/(6l)(sinn6'cos6'Tcosn6lsin6') d9 
■^ Jo 

= - f f{9) sin(n T 1)^ d9 = 26„^i(/), 



1 /■^'^ 

±di(w"'*)-ci(w"'') = - / 2/(6')(±sinn6lsin6' + cosn6'cos6') ^6* 

= - / fie) cos(n T 1)^ d0 = 2a„^i(/). 
I" Jo 



Therefore, we arrive at 



and 



On-U/j = ^^ , On+l(/J = ;^ : 



On-iUJ — :z , ari+il/j — 7i , n > I. 



This simple calculation shows that in order to detect a perturbation that has oscillations 
of order 1/n, one needs to use the first n eigenvectors (e*'^, ^ = 1, . . . , n,) of the Dirichlet-to- 
Neumann operator A/q as Dirichlet boundary data. This is a relatively simple but quite deep 
observation. We conjecture that this result holds for general domains. Another observation 
is that our asymptotic formula is in fact a low-frequency expansion which holds for fixed 
n as e goes to zero. It would be interesting to derive an expansion which is valid for high- 
frequencies, not just for finite n. 



3.3 Proofs of Theorem 13.11 and Lemma 13.31 

We modify uq and ui to the solutions Uq^^ and u^*^ of 



Auf' = 0, 



inR2\B(l-eM,0), 
for 9 e [0, 2n] 
uo{r, 9) — constant + 0(l/r), as r ^ +oo, 



and 



Awf" = 0, 



\S(l-eAf,0), 



uf\l - eM, 9) = -[f{9) + M]drua{l, 9), for 9 G [0, 27r] 



ui{r,9) = constant + 0(l/r), 



where 



M:=max(||/|ji^,||/|U^,l). 
From the fourier expansion of ^, we obtain 



and 

ufir,9)^ 



E 

n=0 



l-eM\» 



an{'^)cosn9 + bn{'9)smn9 



for r > 1 - eM , 



(3.12) 



(3.13) 



-j-oo 



1 -eM\» 



^( —) [an{[f + M]Afo{-^))cosn9 + bn{[f + M]Afo{^))smn9'\. (3.14) 



n=0 



The following is the key lemma to obtain the asymptotic expansion of (u — uq). 



Lemma 3.4 For a 2TT-periodic function ^, we let u and uq he the solution to (|2.5|) and 
p.6p . respectively. 

1. For ^ G C^([0, 27r]), we have the following asymptotic expansion holds uniformly on 
dD,: 

u = u''^^' + euf' + C + 0{€i), (3.15) 

where C is a constant, and O(e^) depends on the Lipschitz constant of f and ||^||,^4. 

2. For a Lipschitz function ^, we have the following asymptotic expansion holds uniformly 
on dDf-. 

O(e^), (3.16) 



,CJ\/ 



where 0{e^) depends on the Lipschitz constant of f and ^. 
Proof. Note that the Dirichlet value of u on 00^ is VP. Using p.l3p and (13.14^ . we obtain 

{u-u^,''~euf')(l + ef,9) 

= C + Y^ \l-( ^^' r Y - ^n{M + f) (a„(*) cosn6l + S„(*) sinnO^ 
■l-eM^ 



n=l 

+ 00 

-E[(^ 

n=l 

C + I + II 



+ e/ 



-1 



a„([/ + Af]AAo(*)) cosnfl + &„([/ + M]AAo(*)) sinvi^) 



where C is a constant. 
Note that 



|1 - (l-t)"| <max{l, 2nt}. 



For *e C4([0,27r]), we have 



|a„(*)|, |S„(^)| < cM^i^ for each n e N, 



(3.17) 
(3.18) 



(3.19) 



and from p.l7p . it follows that 

Now, applying Cauchy-Schwarz inequality, we obtain 



|//|2 < e^ 



1 r/l-eAf\" 



^ n^[\l 



n=l 

+ 00 

X 

ra=l 



ef 



1 



J2 "^^ [""([/ + ^^]-^o(*)) cos 710 + bniif + Af]AAo(*)) sin 710 



<e' 



+c>o 



From dSHD, it follows 

-l-eM\"n2 

n=l 

Therefore we have 



Y- 

L2([0,27r]) ^-^ 11? 
n— 1 



l-eM\" 



1 + e/ 



i:;^['-(t 



+ £,/ 



E;!i(^-)^+Ei^^- 



i<l/e 



i>l/e 



|//M < Ce^ 



(3.20) 



(3.21) 



where C depends on the Lipshitz constant of / and ||^|jc4. 
When ^ is a Lipschtz function, from (|3.20p we have 



(7z-<")(l + e/, 



v-^ r / i — e-M \ "1 / 



n=l 



+ £,/ 

+ 00 



<II^IIl^([o.2.])(E;^[i-(i 



1-eAfN" 



e/ 



2xi 



<Ce5 



D 
Proof of Theorem 13. II For ^ e ^^([0, 27r]), from (|3.15p and the decaying condition of 
M, Ug^^ and m|*^ at infinity, 

u(r, 9) = {u'^"^ + euf ^)(r, 6) + constant + 0(ei), r > 1. 

Let ri be a ball containing D^ , then from the invertibility of the Double layer potential in 
ig(9r2), it follows that 

u{r, 9) = {uf^ + euf''){r, 9) + constant + 0{ei /r), r > 1. 



We calculate that 



+°" (1 - eAf )" - 1 



E 



a„(*) cos nO + &„(*) sin nO 



+00 



^-^ / i — elVl \ " r 



n=0 



[ai(AAo(vI/)/) 



cos6' + 6i(A/'o(^)/)sin6' +C + C>(4 + -), for r > 1, 






depends on the Lipschtz constant of / and 11^1/111^4. 

D 

|a„([/ + M]Uom)\, \bn{[f + M]AAo(*))| < c\, for n G N, 



where C is a constant, and 0{ 
Therefore we prove (13. 4p . 

By the same way, we can prove p.5| 

Proof of Lemma 13^31 Note that 



+00 



XI "''(«»([/ + *^]-^o(*))cosn6' + 6„([/ + M]AAo(*))sinn6') < C, (3.22) 



n=l 



where C depends on jl/jlc^ and jl^jlc*- 
Thus we obtain that 



d_ 



{U-U^0''-<'')\0D. 



+ 00 



E..[i-(t 



l-eAf\» 



+ e.f 



+ 00 



^^iii 



l-eM\" 



n=l 



+ £./ 



+ 00 



e/(^)E"[(t 



- en{M + f) - a„(*)sinn6' + 6„(*)cos7i6' 



a„([/ + M]AAo(*)) sinn0 + S„([/ + Af]AAo(*)) cosn0 



- 1 



l-eAf\" 1 



n=l 

+ OC 



VwEKt 



e/^ 1 + e/ 
l-eA/\" 1 



n=l 



+ e/^ 1 + e/ 



a„('I') cos n6' + 6„(^) sin n6' 
a„([/ + A.f]AAo(*)) cosn^ + 6„([/ + Af]AAo(*)) sinn^ 



0(6^). 



(3.23) 



There exists a constant C which depends on the Lipschitz character of dD^, see [T], such 



that 



d_ 
dv 



{u-uf' -euf') 



L2(a_D^ 



<C 



d_ 
df 



(M-<^^-eufO 



L^{dD,) 



where T is the unit tangent vector on dD^. From (|3.23p and the fact that -gg = [1 + ef - 
0(6^))^, it follows that 



|-(.-<--..r; 



LHdD,) 



0{e-^). 



(3.24) 



From dnH), (EZD and ([XM]) . we have 



= (1 - e/)(VK^^ + 6<0,iVe) + 0{ei) 



d_ 
dr 

^( 
i9r 



("r+e<o 



r=l+e/(e) (90 ^ 



r=l+e/(e) 



+ 0(e^) 



£j\/ 



eu 



;Af^ 



r=l+e/(e) 



e/^ + 0(e2). 



(3.25) 



We compute 

+ 00 



r=l+e/(e) 



^ — — jf-TT^r) (-'^)(an(*)sinn0 + 6„(*)cosn6' 



e/ V 1 + e/ 



+ 00 



Y. Y^f {tTTtY''^^''^^^ + MWoW) cosnO + fe„([/ + MWoi^)) sinnO 



Since 



.— — (— — r) =-n + en(n+l)/ + en2M + 0(eV), 
1 + £/ V 1 + e / y 



1 /l-eM\" 
^nz 1 -^ r 1 = en + £71 



1 + e/ V 1 + e/ 



l-eM\" 



1 + e/ 



-1 



+ 0(e2n), 



we have 



li<' + ^0 



r=l+ef{e) 



M)(*) + e(l?oM)(*)/ - AAo(AAo(vl/)/)) + 0{e 



From p.25p . we prove the lemma. 



n 



4 Formulation of the Acoustic Problem 

Analogously to the Laplacian one, we study the inverse scattering problem of reconstructing 
a sound-soft obstacle, call it D^, whose boundary is the perturbation of the unit circle and 
is given as (12. ip . 



4.1 Inverse Scattering Problem 

For a incident field v\ we denote v'^ and Vq as the scattered field from D^ and D, respectively. 
In other words, v'^ and Vq are the solutions to 



r Av' + k^v' = 0, 



d_ 
dr 



on dDf, 



v^{r,9) — ikv'^{r,9) — o{r ^), r — > +c 



(4.1) 
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and 

Aug + k^v'o = 0, in M2 \ jj^ 

Vq = ~v\ on dD, ^42) 

—v^{r,e)~ikv^{r,e)^o{r-^), r — » +00, 

Here we used the polar coordinates a;(r, 9) = (r cos 9, r sin0), and the wave number k is given 
by a positive constant. 

4.2 Fixed Dirichlet boundary value problem 

For a 27r-periodic continuous function ^P, we let u be the solution to the Helmholtz problem 
with the prescribed boundary data ^ on dD^, i.e., 

Au + k'^u = 0, mR'^\D,, 

u{l + ef{9),9) = ^(9), for 9 G [0, 27r], (43) 

-—u(r,9)^iku(r,6) = o{r^^), r — > +00. 
or 

The solution u corresponding to the unit disk D satisfies that 

Auo + fc^uo = 0, in M^ \ D, 

Mo(l, 6*) = *(6l), for^e [0,27r], 



—uo{r,9)~ikuo{r,9)~o{r^'^), r — > +00. 



(4.4) 



We investigate the Far-Field difference between u and uq, especially when ^ is a C^- 
function or a Lipschitz function. 

5 Acoustic far-field formula and inversion algorithm 

5.1 Asymptotic Far-field expansion for the Dirichlet problem 

We parametrize the unit circle dD by 9 E [0, 27r] and expand ^ as 

*(0) = ^c„(*)e^"^ 

nez 

where c„(^) is the fourier coefficient with respect to e™^. By the uniqueness of the exterior 
Dirichlet problem, it follows that 






"«(^'^) = E TTlTTir ^"(*)^"'- (5-1) 



Define the Dirichlet-to-Neumann operator JVq with respect to D by 

A/'o : uo\s -^ drUols, 
then, in a pseudodifferential fashion, Ao can be written as follows (see [7]): 

Afoi^m = ^al(n,fc)c„(*)e™^ (5.2) 



11 



where the so-called discrete symbol cti is given by 

(Ti{n, k) = k — "^ = -k — ^^- h \n\. 

Thus, for fixed k, we have 

ai{n,k) ~ |n|, as \n\ -^ oo. (5-3) 

By the same way as the electric problem, we can consider ui as the solution to (|4.4|) with 
the boundary value (— A/'o(^)/) on dD instead of ^, and it follows 

Hl'hkr) 






ui(r, U} = - } ^ . (1),.. c„^7Vo(,WJ/je 



It is known that, for a fixed n, the Hankel function of the first kind satisfies 



-^H (2^) ^ vJ^'^""^~'"'^^ + 0(|x|-i), X » 



(5.4) 



We refer to [3] for more properties of the Hankel function. 
Choose iV G N satisfying that 



Y.\cn{fm^))\=0{e^), (5.5) 

\n\>N 



then we have the following lemma. More precise proof is given in the Subsection 15.31 
Theorem 5.1 1. Let ^ G C'"'([0, 2ti]) and u he the solution to (|4.3[) . For r ^ I, we have 

{u-uo){r,9) = -e\[^e'>'^ V ^"(•^o('^).^) e-»(f + ^)e'"^+0(el/V7+e/r). (5.6) 

\n\<N "-\n\ y^' 

where N is defined by (|5.5p . and 0(e^/-\/r + e/r) depends on the Lipschitz constant of 
f and||*||c4. 



^. Let ^ 6e a Lipschitz function and u be the solution to (|4.3p . FFe have that 

{u-uo){r,9)^0{e^/V^), /orr»l, (5.7) 

where 0(t^ j \fr) depends on the Lipschitz constant of f and ^. 

We define 

Voi^m := ^ a^in, fc)c„(vI/)e™^ (5.8) 

nez 
with 

Lemma 5.2 For / e C2([0,27r]) and * G C"*([0, 27r]), we have that 

Kf^ = AAo(*) + eAA/(*) + 0(ei), 
w/iere 
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5.2 Algorithm for the Inverse Shape Problem 

Let w' be the incoming wave, and define v'^ and Vq as the solution to (|4.ip and (|4.2p . 
respectively. Note that 

v\l + ef{0),0) = vHl, 9) + ef{B)drv\l, 0) + 0{e^). 

Applying Theorem 15.11 by letting 'f = —v'^{l,d) and ^ = —ef{6)drV^{l,9), we have for 
r > 1 that 



g~i(f + ^)gme 



fY ^^ CnifAfoiv'lao)- fdrV'loD 



r(l) 



n|<A' 



Hr {k) 



where N is defined by X)|n|>Af c,i( /A/'o(u*|a_D) — /Sri^'laDJ = 0{e-^). This yields to stable 

reconstruction of the Fourier coefficients c„(/A/'o(w*|aD) — /9rW*|aD 1 for n such that H^!{k) 
is not too big. 

Suppose now that v'^ satisfies 

then by measuring ci(v''^ ~ ^o)' then we can reconstruct Cm{f). 



5.3 Proofs of Theorem 15.11 and Lemma 15.21 

We modify mq and ui as Uq^' and w|*^ which satisfy 

(A + k'^)uf' = 0, in K^ \ B(l-eM,0), 

w5*^(l - eM, 6*) = ^(61), for 9 G [0, 27r] 

^ul^\r,9) - ikul^\r,9) = o(r-3), r > +oo. 



and 



f (A + fc^)Mf ^ = 



inM^\S(l-eAf,0), 



uf^(l - eM, 9) = -[f{9) + M]AAo(*), for 9 e [0, 27r], 



or 



-oo, 



where M is the constant defined by p.l2p . From the fourier expansion of ^ and the unique- 
ness of the exterior Dirichlet problem, Uq^^ and uf*^ have the expansion as follows: 



r(i) 



^o"(^^) = Et7(T 



Hnikr) 



t:^^H\^>{k-ekM) 



c„(*)e™^ 



ur(r,« 



^_ H[^]{kr) 



£„([/ + M]AAo(*))e' 



'in6' 



«a^H(fc-efcAf) 
We have the following key lemma to prove Theorem 15.11 and Lemma 



(5.9) 



(5.10) 
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Lemma 5.3 1. For '^ G C^([0, 27r]), we have the following asymptotic expansion for the 
solution u to (|4.3p holds uniformly on dD^: 



where 0{e^) depends on the Lipschitz constant of f and W^^Wc^- 



(5.11) 



2. For a Lipschitz function "^jWe have the following asymptotic expansion for the solution 
u to (14. 3p holds uniformly on dD^: 



where 0{e^) depends on the Lipschitz constant of f and ^. 
Proof. From (|5.9p . (|5.10p and the boundary condition of u on dD^, we have that 



(5.12) 



^^Ehi 



nez ^^|„| 



Hi:^>{k + ekf) 
^HJ'Uk-ekM) 



1 



c„([/ + M]AAo(*))e™^ 



and 



^_r HJ-^hk + ekf) 



n& 



H^ik-ekM) 



c„(*)e' 



Note that 



<(fc + ri)-<(fc) 



<eiii m, 



1)' 



|„| IU~([fe-e|t|, fc+e|t|]), 



rWn. , .^^ it(i)^m ..E7(l)V;.^ ^ ^ ^ iiu(i)"i 



From the fact that 



.^2 

2 



|n| V / l^+l| ^ ^ y |n| ^ /' 



we can show that 



ijflVfc + efc/) 

1 7TT = 0(en), 

ijf„^l'(fc-efcM) 

H^^^ [k + efcf) 

iJ«(fc~6A;M) 



(5.13) 
(5.14) 



where 0(en) and O(e^n^) depend on Af and k. Moreover, \Il\}\[z') is decreasing function 



for z > 0, and 






o(i). 



(5.15) 



Using (I5.13p . (|5.14p and (I5.15p . we can prove the lemma by the same way to prove Lemma 

D 
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Proof of Theorem 15. II At first, we assume ^ e C*. From the solution expression using 
boundary integral methods (for example, see [5j), we can show that 



iu-u^^^' - e<")(r,0) = ^h- <" - euf'WmoD,), as r 



-oo. 



(5.16) 



Therefore 



{u - wo)(r, e) ^ (w^*'^ + €uf' - u^){r, 6) + 0(e*/V^), as r ^ +c 



From ([5?T|) . ([O]) and (jS.lOp . we compute that 






Hi:.iikr) r H\:^>{k) 



^fi-«(fc) LFW(fc-efcM) 



r(i) 



N 

rr(l 



4^1 (fc'^) 



(1)/ 



ne 



tHi;;{k)^ Hi;>{k~ekM) 



2^ rr(l) 






H 

rr(l 



(1)/ 



^m:(fc)L m;(fc-eMf) 



<?(^) 



■I + 11 + III. 



fi„(^)e"' 

Cn(/AAo(*)) + c„(AfAAo(*)))e''"« 
-l-eMcri(n,fc)j c„(^)e"^ 
c„(/AAo(vI/))e™^ 
Mc„(AAo(vI/))e™^ 



From (|5.13p and (|5.14p with replacing / by and (|5.16p . it follows 



I + 111= -^0{e^). 



Using (|5.4p and (|5.16p . we obtain 



//= -€\l—e 
■nr 



2 ... ^ cMKm^-.,.^l^)^.^e^lo{,)+0{eilV-r). 



ikr 



E 

n|<Af 



<?(fc) 



(5.17) 



(5.18) 



When $ is a Lipschitz function, we have 

{u-uo){r,0) = (u^*^-uo)(r-,6')+0(e5/V^) = 0(eVV^), as r ^ +oo. 



n 



Proof of Lemma 15.21 Note that dD is a C^-domain, and using boundary integral 
methods, we have that (see [5]) 






^^' euf'] 



cO'°(an,) 



<C 



u~ u, 



eM ^.M 



ci.°(an,) 



(5.19) 



By the same way as the conductivity case, we can prove the lemma by calculating ^("o^^ - 



eu\^) 



D 
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